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SCOPE 


This document describes the acquisition of daily weather data from North Dakota and 
Kansas,, the reduction of U.S.S.R. synoptic observations (4-8 observations daily) to 
dally data and conversion of all three data sets to the formats utilized by ad- 
vanced wheat yield models,. Data from the three areas are available on tape, those 
of each area having a unique format. Tapes were built by converting these data 
to one of two possible formats, depending on the amount of data available. 

North Dakota and Kansas data contain daily temperature maximums and mlnimums, as 
well as total daily precipitation. North Dakota precipitation Information is 
less detailed than that of Kansas. No significant manipulation is required to 
convert these data to the needed formats. 

Synoptic data are available for the U.S.S.R., with a maximum of eight 3-hourly 
observations per day containing detailed weather information recorded at the 
time of observation. An algorithm is applied to the U.S.S.R. data to estimate 
dally maximum and minimum temperatures and sum the daily precipitation. 

Nineteen tapes have been obtained to date from NOAA's National Climatic Center 
(NCC) at Asheville, North Carolina including synoptic weather data from WMO blocks 
33, 28, and 35. In general Block 33 covers the western three-fourths of the 
Ukraine, Block 28 covers the area on both sides of the Ural Mountains and Block 35 
includes the area between the Ural Mountains and the Aral Sea. Within the wheat 
areas in these blocks daily precipitation and max-min temperature records have 
been derived for forty-seven stations for most years during the period 1965 
through 1975. Additional tapes are available from NCC from which similar records 
could be derived for the roughly RO percent of the U.S.S.R. not included in 
Blocks 33, 28, and 35. 


2. APPLICABLE DOCUMENTS 


• Job Order 63-1347—4963 

* Action Documentation 63-1347-4963-08 (Task II, AD if 1) 


3. SYSTEM DESCRIPTION 


3 « i HARDWARE DESCRIPTION 

Those data conversion programs will operate on an IBM 360/370 system with 
Fortran and PL/I compilers and magnetic tape drive. In the event that a 
PL/I compiler is not available, one short program can be rewritten in Fortran. 
These programs can also operate on a PDP 11/45 batch environment with a 
Fortran compiler and magnetic tape drive. 


3 • 3 NORTH DAKOTA DATA BASE FORMAT AND CONVERSION PROGRAM 

The tape obtained from North Dakota State University contains dally data for 
59 North Dakota weather stations, each station on a separate file. HDSU 
personnel had assembled these data in 26-character records, with each record 
containing weather data for one day/station. 

The North Dakota 26-character format is different from the standard 26-character 
format used for the advanced wheat yield models (Format I, referred to by Kansas 
State University personnel as "modified Kansas format") , A short PL/I program 
was written to convert the data of 12 stations to Format I. This program can 
be rewritten In Fortran if a PL/I compiler is not available. Documentation is 
given in Appendix A. 


3 . 3 KANSAS DATA BASE FORMAT 


The tape obtained from Kansas State University contains daily data for 39 Kansas 
weather stations, with each station on a separate file, KSU, in turn, had 
obtained their data from the National Climatic Center in Asheville, North 
Carolina, and had filled in missing data from local sources. 

The tape format is referred to by KSU personnel as "Format II," It consists 
of 960-character records, with each record containing the necessary input 
weather data for one month/station. Documentation is given in Appendix B. 


3.4 U.S.S.R. DATA BASE FORMAT AND CONVERSION PROGRAM 

Nineteen tapes were obtained from the National Climatic Center in Asheville, 
North Carolina, covering World Meteorological Organization (WMO) blocks 28, 

33, and 35 for the years 1965-1975, These tapes contain synoptic data for all 
reporting weather stations within these three areas, recorded in WMO tape deck 
format 9605 . 

The program to perform this conversion (max-min algorithm) and its subroutines 
are described subsequently. 


3.4.1 MAX-MIN ALGORITHM 


The Statistics Department of Kansas State University has developed a non-linear 
regression model to estimate maximum and minimum temperatures from synoptic data. 
This algorithm was programmed in Fortran by LEG personnel, using a main program 
supported by seventeen different subroutines and functional subprograms. The 
program is unique, with no other major library subroutines or linkage needed. 

Each of the subroutines will be described below under separate paragraph 
headings; flowcharts and sample listings of the main program and all subroutines 
will be found in Appendix C. The results are stored on a tape in Format I, 

Program MINIMAX is designed to estimate the daily maximum and minimum temperatures 
for weather stations in the U.S.S.R. with greater than 90% accuracy when compared 
with Soviet-published daily weather data. 

The algorithm is based on temperature readings recorded every three hours over a 
24-hour interval. Two possible situations are associated with the given 
observations; 

1. Fewer than eight readings in a 24-hour interval. 

The estimation for minimum and maximum temperatures will be based 
on specified readings. Suitable correction factors derived from 
the algorithm are applied to the selected temperature readings to 
determine the day’s minimum and maximum temperatures, 


The daily 

minimum will 

utilize the following 

temperature 

readings: 





Time 


Hours 



Zone 

1st Choice 

2nd Choice 

3rd 

Choice 

" 00 

03 

06 


00 

01 

03 

00 


06 

02 

03 

00 


06 

The daily 

maximum will 

use the following 

temperature readings: 

Time 


Hours 



Zone 

1st Choice 

2nd Choice 

3rd 

Choice 

00 

15 

12 


09 

01 

12 

15 


09 

02 

12 

09 


15 


2. Eight readings in a 24 -hour interval. 

The highest and lowest readings among the eight, with proper correction 
factors, will be the daily maximum and minimum temperatures, respectively. 


3.4.2 SUBROUTINE CONV 

The subroutine CONV is used to convert special characters (missing data) from 
the temperature and precipitation variables to numeric values. Input and output 
arguments are TEMP and PRECIP, 
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3.4.3 SUBROUTINE CMDTOJ 


This subroutine converts the input arguments YR, MONTH and DAY into the Julian 
date and output argument, JILIAN. 


3.4.4 SUBROUTINE EST1 

If fewer than eight observations for a particular day have been recorded, this 
subroutine chooses the best possible hour and temperature reading to use in 
calculation of the sun time (ST) for the day's minimum temperature. If all 
eight observations have been recorded, this subroutine is not utilized. The 
input arguments are the station's time zone HO, an array of the day's temperature 
observations KHRTEM, the sun time correction factor SUNCOR, and the hour of 
observation HOUR. Output arguments are the sun time ST, and the position of 
the chosen temperature reading within the array of a day's observations HT. 


3.4.5 SUBROUTINE EST2 

If fewer than eight observations for a day have been recorded, this subroutine 
chooses the best possible hour and temperature reading to uBe in calculation of 
the sun time (ST) for the day's maximum temperature. If all eight observations 
have been recorded, this subroutine is not used. Input and output arguments are 
identical to those used in ESTL. 


3.4.6 SUBROUTINE ESTDTN 

If fewer than eight observations have been recorded, this subroutine estimates 
the value for DTW, the correction factor for the daily minimum temperature. 
Regression coefficients were provided by che KSU Statistics Department. Input 
arguments are an array of regression coefficients B, an array of variable 
combinations X, the sun time ST, daylength DL and the average monthly temperature 
range for the station TR; the output argument is DTN. 


3.4.7 SUBROUTINE ESTDTX 

If fewer than eight observations have been recorded, this subroutine estimates 
the value for DTX, the correction factor for the daily maximum temperature. 
Regression coefficients were provided by the KSU Statistics Department. Input 
arguments are the same as those used in ESTDTN; the output argument is DTX. 


3.4.8 SUBROUTINE MAXMIN 

If fewer than eight observations are recorded, this subroutine estimates the 
daily maximum and daily minimum temperatures from the specified hour temperature 
readings. Input arguments are the correction factor for the daily maximum 
temperature DTX, the correction factor for the minimum temperature DTN, the 
recorded temperature at the hour chosen to calculate the maximum HTDX, and the 
hour chosen to calculate the minimum HTDN. Output arguments are the day's 
estimated maximum and minimum temperatures, DMAX and DMIN, respectively. 
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3.4.9 SUBROUTINE MAX 


When all eight temperature observations have been recorded, this subroutine Is 
used to find the highest temperature of the eight. The Input arguments are an 
array of the day's observations KURTEM, an array of subscripts to KHRTEM called 
K, and the sun correction factor SUNCQR. Output arguments are the maximum 
observed temperature OBSMAX and the hour at which this max occurred b*fiX. 


3.4.10 SUBROUTINE MIN 

When all eight temperature observations have been recorded, this subroutine Is 
used to find the lowest temperature of the eight. The Input arguments are the 
same as those for subroutine MAX; the output arguments are the minimum observed 
temperature OBSMIN and the hour at which this minimum occurred STMN. 


3.4.11 FUNCTION DLNGTH 

This functional subprogram estimates the daylength of each Julian day. The 
input arguments are the latitude of the station being considered XLAT and the 
Julian date DATE. The output is transferred bock to the main program by the 
name of the funcltonal subprogram DLNGTH. 


3,4.12 SUBROUTINE PSUD01 

When all eight temperature observations have been recorded, this subroutine 
estimates the correction factor for the daily minimum temperature DTN for 
stations located in time zone 00. Input arguments are an array of regression 
coefficients B, an array of variable combinations X, the hour at which the 
minimum temperature occurred STMN, daylength DL, average monthly temperature 
range TR and the sun time correction factor SC; the output argument is DTN. 


3.4.13 SUBROUTINE PSUD02 

When all eight temperature observations have been recorded, this subroutine 
estimates the correction factor for the daily maximum temperature DTX for 
stations located in time zone 00. Input arguments are an array of regression 
coefficients B, an array of variable combinations X, the hour at which the 
maximum temperature occurred STMX, daylength DL, average monthly temperature 
range TR and the sun time correction factor SC; the output argument is DTX. 


3.4.14 SUBROUTINE PSUD03 

When all eight temperature observations have been recorded, this subroutine 
estimates the correction factor for the daily minimum temperature DTN for 
stations located in time zone 01. Input and output arguments are identical 
to those of PSUD01. 
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3.4.15 SUBROUTINE PSUD04 


When all eight temperature observations have been rticftf-jerf, this subroutine 
estimates the DTX for stations in time zone 01. Input and output arguments 
are identical to those of PSUD02. 


3.4.16 SUBROUTINE PSUD05 

When all eight temperature observations have been recorded, this subroutine 
estimates the DTN for stations in time zone 02, Input and output arguments 
are the same as those of PSUC01. 


3,4.17 SUBROUTINE PSUD06 

When all eight temperature observations have been recorded, this subroutine 
estimates the DTX for stations in time zone 02, Input and output arguments 
are the same as those of PSUD02, 


3.4.18 SUBROUTINE OUT 

This subroutine is used to put tiw output information into Format I (see 
Appendix A) . The input arguments are the station number STATN, year SAVYR, 
month SAVMON, day SAVDAY, estimated maximum temperature DMAX, estimated 
minimum temperature DMIN, the day's precipitation total SUMPRE . Output 
arguments are STATN, SAVYR, SAVMON, SAVDAY, maximum temperature MAX, 
minimum temperature MIN and the day's precipitation PRE, 
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4. OPERATION 


To operate program MIftTMAX the number of stations to be used, NUMSTA, must be 
declared. The following card input is necessary, using one card/station 5 


Columns 

6 

7 - 10 

11 - 70 

11-15 

16-20 


Variable name, description 
STNARR, station number 
blank 

SUNARR, sun correction factor in thousandths of degrees 
longitude, no decimal point 

TRARR, monthly temperature ranges, in hundredths °F, 
no decimal points 

January temperature range 
February temperature range 


66-70 December temperature range 

71 blank 

72 - 75 LATARR, station latitude in hundredths of degrees, 

no decimal point 

76 - 77 TZARR, station time zone location 

78 - 80 blank 

If more than 30 stations will be used, space for the above arrays must be 
increased. Output data from MINIMAX will be in Format I. 


After conversion to either Format I or Format II, the North Dakota, Kansas and 
U.S.S.R. data bases may be directly accessed by the advanced wheat yield models. 
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APPENDIX A 
FORMAT I 


FORMAT I 


Character String 
(per record) 

1 - 6 
7-8 
9-10 
11 - 12 

13 

14 - 16 
17 - 19 
20 - 22 
23 - 26 


Station number 

Year 

Month 

Day 

Blank 

Max Temperature * 
Min Temperature * 
Blank 

Precipitation ** 


* Temperatures are in whole degrees Fahrenheit with a leading minus sign if 
negative; missing temperatures are indicated by blanks or 999. 

** Precipitation is in hundredths of inches with no decimal point; missing 
precipitation is indicated by 9999. 
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APPENDIX B 


FORMAT II 


FORMAT II 


Character String 
(per record) 

1 - 2 
3-6 
7 - 8 
9-10 
11 - 134 
135 ~ 227 
228 - 320 
321 - 382 
383 - 413 
414 - 475 
476 - 568 
569 - 661 
662 - 785 
786 - 878 
879 - 940 
941 - 960 


State Index // 

Station Indes 

Year 

Month 

Precipitation (31 values) 

Max Temperature ,13 values) 

Min Temperature (U values) 

Max water temperature - Pan (31 values) 
Estimated precipitation (31 values) 

Snowfall (31 values) 

Snow Depth on ground (31 values) 

Water equivalent of snow on ground (31 values) 
Wind movement (31 values) 

Evaporation (31 values) 

Min water temperature - Pan (31 values) 

All l's 


Precipitation (fields of 4) 

Characters 11 - 134, 31 values, field of 4. For months with less than 31 days 
fill remaining fields with 9999. 

Contents of each field 


0001 - 9997 00.01 - 99.97 inches to hundredths 

0000 for trace 

0000 for no precipitation 

9998 amount included in subsequent measurement 

9999 missing or not reported 

Maxi mum Te mperature (fields of 3) 

Characters 135 - 227, 31 values, field of 3. For months with less than 31 days 
fill remaining fielus with 999. 

Contents of each field 

001 - 299 -99. to 199. °F. Bias each value by +100°F thus 

character string will contain values from 1 to 299 

999 not reported or missing 

Minimum Temperature (fields of 3) 

Characters 228 - 320 

Contents of each field (Same as for Maximum Temperatures) 
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Maximum Water Temperature - Pan (fields of 2) 

Characters 321 - 382, 31 values, field of 2. For months with less than 31 days, 
remaining fields are 99, 

Contents of each field 

00 - 98 30 to 128°F Each temperature value is biased by -30°F; if tempera- 

ture is greater than 128°F it is set to 128°F, if 
temperature is less than 30°F it is set to 30°F. 

99 Missing or not reported 

Estimated Precipitation (fields of 1) 

Characters 383 - 413, 31 values, field of 1. For months with less than 31 days, 
remaining fields are 9, 

Contents of each field 

9 When daily precipitation was not estimated or when it 

is missing or nor reported. 

i. When daily precipitation was estimated. 

Snowfall (fields of 2) 

Characters 414 - 475, 31 values, field of 2. For months with less than 31 days, 
fill remaining fields with 99. 

Contents of each field 

01-97 1 to 97 Inches rounded to nearest whole inch (if over 

97 inches change to 97) 

00 No snowfall 

98 Amount included in subsequent measurement 

00 Trace 

99 Missing or not reported 

Snow Depth on Ground (fields of 3) 

Characters 476 - 568, 31 values, field of 3. For months with less than 31 days, 
fill remaining fields with 999. 

Contents of each field 


001 - 998 

000 

000 

999 


1-998 Inches (whole inches) 

None 

Trace 

Missing or not reported 


Water Equivalent of Snow on Ground (fields of 3) 

Characters 569 - 661, 31 values, field of 3. For months with less than 31 days, 
fill remaining fields with 999. 

Contents of each field 


000 - 998 
999 


00.0 to 99.8 inches of water 
Missing or not reported 


*N0TE. .. .water equivalent data available only at first-order stations beginning 
October, 1963. 
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Wind Movement (fields of 4) 

Characters 662 - 785, 31 values, field of 4. For months with less than 31 days, 
fill remaining fields with 9999. 

Contents of each field 

0000 - 9997 Whole miles (if over 9997 miles, change to 9997). 

9998 Amount included in subsequent measurement 

9999 Missing or not reported 

Evaporation (fields of 3) 

Characters 786 - 878, 31 values, field of 3, For months with less than 31 days, 
fill remaining fields with 999. 

Contents of each field 

001 - 997 0.01 - 9.97 (to hundredths) 

998 Amount included in subsequent measurement 

999 Missing or not reported 

Minimum Water Temperature - Pan (fields of 2) 

Characters 879 - 940 

Contents of each field (same as for Maximum Water Temperature - Pan) . 
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